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ABSTRACT 

Based  on  a  theoretical  analysis  of  the  cable  and  payload  dynamics  during 
lowering  or  raising  heavy  loads  in  the  deep  ocean  given  in  Project  Trident  Technical 
Report  No.  1370863,  further  calculations  of  the  maximum  dynamic  stresses  expected 
in  the  lowering  cable  are  presented  covering  a  wide  range  of  cable  and  payload 
parameters.  The  theoretical  analysis  is  adapted  to  a  proposed  design  procedure,  and 
two  typical  design  examples  are  given,  the  results  of  which  are  discussed  in  terms  of 
the  safety  of  the  lowering  or  raising  operations. 

In  order  to  make  the  design  procedure  applicable  with  a  greater  degree  of 
confidence,  it  is  considered  necessary  to  make  measurements  of  cable  tensions  and 
load  and  ship  motions  during  a  full-scale  operation  to  fill  in  deficiencies  of  data 
and  provide  a  basis  for  verification  of  theory  and  calculations.  In  particular,  data 
are  needed  on  the  coefficients  of  drag  and  mass,  which  at  this  stage  must  of  necessity 
be  estimates. 
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A 


INTRODUCTION 


The  work  described  herein  was  carried  out  os  part  of  BuDocks  Task  No. 
Y-F015-01-01-001,  "Structures  in  Deep  Ocean,"  which  originated  from  the  require¬ 
ment  of  the  Bureau  to  attain  a  deep  ocean  engineering  capability  in  keeping  with  the 
increased  emphasis  on  the  deep  ocean  as  an  operating  environment  for  naval  forces. 
This  report  is  the  result  of  work  performed  under  Task  No.  Y-F015-01  -01  -001(b), 
"Mechanics  of  Raising  and  Lowering  Heavy  Loads  in  the  Deep  Ocean." 

The  objectives  of  this  work  were  to  analyze  and  report  on  the  results  of 
predictions  of  the  forces  in  lines  and  the  acceleration  and  displacements  of  loads  of 
various  shapes  during  raising  and  lowering  operations  in  the  deep  ocean.  The  study 
is  based  on,  and  an  extension  of,  a  theoretical  analysis  by  Arthur  D.  Little,  Inc., 
given  in  Project  Trident  Technical  Report  No.  1370863  of  the  Bureau  of  Ships, 
entitled  "Stress  Analysis  of  Ship-Suspended  Heavily  Loaded  Cables  for  Deep  Under¬ 
water  Emplacements. "  ^ 


THEORY 


The  problem  considered  in  this  report  is  that  of  a  load  suspended  from  a  ship  or 
moored  platform  by  means  of  a  single  cable,  as  shown  in  Figure  1.*  The  maximum  depth 
for  such  a  lowering  or  raising  operation  is  assumed  to  be  on  the  order  of  20,  000  feet. 

The  requirements  are  (1)  to  analyze  this  problem  so  as  to  predict  the  cable  and 
the  load  dynamics,  in  particular  the  maximum  dynamic  stress  induced  in  the  cable  as 
a  result  of  the  motions  of  the  suspension  point,  and  (2)  to  provide  a  design  procedure 
for  evaluating  such  stresses  for  a  given  load  under  specified  conditions  of  sea-surface 
oscillations. 

Solution  of  this  problem  is  recognized  to  be  difficult  in  view  of  the  nonlinearities 
introduced  in  the  damping  due  to  drag  forces  of  vertical  oscillations  of  the  load.  A 
simplified  solution  has  been  obtained  by  Arthur  D.  Little,  Inc.,  as  given  in  the  Project 
Trident  report.  A  brief  resume  of  the  analysis  presented  in  that  report  is  given  here, 
details  of  which  may  be  found  in  Appendix  A. 

The  equation  or  morion  of  an  element  of  the  cable  initially  located  at  a  distance 
x  from  the  support  point  is  given  by 
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*  Figures  1  through  35  are  presented  immediately  after  the  main  body  of  text. 


The  boundary  condit'rm  at  the  load  is  given  in  the  form 
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A  secondary  boundary  condition  is  that  the  displacement  of  the  suspension  point  is 
known  for  all  time. 

In  the  above  equations,  u(x,t)  is  the  displacement  at  time  t  of  a  point  on  the 
cable  originally  located  at  a  distance  x  from  the  suspension  point.  Hence  the  second 
boundary  condition  is  that 


u(o,  t)  is  known  for  t  >  0. 
The  following  nomenclature  applies: 
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where*  p 

c 

s 

E 
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density  of  the  cable 

material  cross-sectional  area  of  the  cable 
modulus  of  elasticity  of  the  cable 

constant  of  friction  on  the  cable  due  to  the  surrounding  water 


*  Notations  are  defined  where  they  first  appear  and  are  summarized  for  convenience 
on  a  foldout  page  at  the  back  of  the  report. 
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c  =  velocity  of  sound  in  the  cable 
MQ  =  virtual  mass  of  the  load 
Cq  =  drag  coefficient  appropriate  to  the  load 
A  =  projected  area  of  the  load  in  the  direction  of  motion 
L  =  length  of  cable 
p  =  density  of  sea  water 

The  difficulty  in  obtaining  an  exact  solution  to  Equation  1  subject  to  the 
boundary  conditions.  Equations  2  and  3,  arises  from  the  nonlinear  term  | du/ot*  | (bu/St1 ) 
in  Equation  2„  This  difficulty  is  avoided  in  the  A.  D.  Little  report^  by  an  approxi¬ 
mation  which  is  described  in  Appendix  A. 

Defining  a  normalized  displacement  amplitude  U'  equal  to  U  divided  by  J  U0  j 
and  noting  that  Uj  is  the  value  of  U1  at  the  load,  a  solution  for  U'  as  a  function  of 
x'  is  given  by 


U'  -  Ujcostu'y'  +  Csinou'y'  (5) 

where  C  is  a  complex  constant  and 


y'  =  1  -  x' 

(6a) 

,  *L 

uo'  =  — 

c 

(6b) 

Hence  the  maximum  value  of  the  dynamic  stress  in  the  cable  is  given  by 


(£] _ =  (u;1)^  (U  j)^  [  1  +  tan(p(tan^  +  sec'I') 
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In  the  A.  D.  Little  report,  Equation  7  together  with  Equations  8,  9,  and  10 
was  solved  by  use  of  a  digital  computer  to  give  the  maximum  dynamic  stress  as  a 
function  of  the  nondimensional  frequency  co'  for  various  values  of  /3  and  4.  A 
similar  procedure  was  adopted  ir  this  report  for  two  reasons,  firstly  to  investigate 
the  variation  of  stress  over  a  wide  range  of  a;',  0,  and  p,  and  secondly  as  a  means 
of  providing  the  basis  of  a  design  procedure  for  cables  used  for  lowering  or  raising 
heavy  loads  to  or  from  the  deep  ocean  floor. 

The  cable  and  load  system  considered  herein  is  a  part  of  the  overall  lowering 
system  consisting  of  the  vessel  from  which  the  operation  is  performed,  its  response 
to  the  wave  action  present  during  the  lowering  or  raising  process,  and  the  resulting 
oscillations  of  the  cable  and  load.  Within  existing  theoretical  limitations  of  knowl¬ 
edge  about  waves  and  ship  motions,  and  under  the  restrictions  of  a  linear  theory,  the 
problem  of  the  response  of  a  ship  or  platform  to  a  particular  sea  state  has  been  solved 
in  terms  of  certain  probabilistic  models  by  Kaplan  and  Putz.^  Pierson  and  Holmes^ 
in  a  note  on  the  engineering  applications  of  the  Kaplan  and  Putz  report  outlined  a 
procedure  for  the  determination  of  the  response  of  a  drilling  barge  to  sea  states  3,  4, 
and  5.  The  results  are  obtained  in  terms  of  the  probability  of  occurrence  of  various 
amplitudes  of  motion  in  heave,  surge,  sway,  yaw,  pitch,  and  roll.  The  Cuss-I 
ocean-bottom  drilling  barge  was  used  as  an  example,  but  the  calculations  os  carried 
out  by  Kaplan  and  Putz  may  be  applied  to  other  ships  or  moored  platforms,  given  the 
use  of  a  digital  computer. 

Details  of  the  program,  which  was  written  for  an  IBM  1620  computer,  are 
given  in  Appendix  C.  For  the  purposes  of  this  analysis,  the  calculations  were 
divided  into  sections  based  on  the  relative  values  of  u;'  required  for  prototype 
computations. 

Equation  6b  relates  the  required  range  of  a;'  to  the  length  of  and  velocity  of 
sound  in  the  cable.  It  is  assumed  that  x,  the  frequency  of  oscillation  of  the  cable- 
suspension  point,  has  a  maximum  value  on  the  order  of  2. 00  radians  per  second  and 
that  the  maximum  length  of  the  cable  is  20,000  feet.  Then  the  required  range  of  w* 
is  determined  by  the  velocity  of  sound,  c,  in  the  cable.  For  steel  and  polypropylene 
cables,  c  is  approximately  12,000  and  2,  000  feet  per  second  respectively,  resulting 
in  maximum  nondimensional  frequencies  of  3.33  and  20.00. 
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Initially  the  range  of  $  was  chosen  as  0. 10  to  7,  00,  and  values  of  u  equal  to 
0.10,  0.50,  1.00,  2.00,  5.00,  and  10. 00  were  used.  Preliminary  computations 
for  a  typical  design  problem  indicated  that  lower  values  of  p  would  also  be  required, 
and  corresponding  additional  computations  were  carried  out  as  shown  in  the  results 
which  follow. 

Pierson  and  Holmes  indicate  the  methods  whereby  the  root-mean-square  (R./S) 
values  of  motion  in  each  mode  for  each  sea  state  may  be  determined.  For  the  purposes 
of  this  report,  the  oscillation  of  most  concern  is  that  in  heave,  and  knowing  the  RMS 
value  of  heave  motions  in  sea  state  4,  for  example,  estimates  can  be  made  of  the 
extreme  value  of  heave  to  be  expected  in  a  given  time.  This  procedure  thus  provides 
a  basis  for  specifying  the  range  of  |  U0|  to  be  used  in  determining  ,'c!  jes  for  use  in 
the  design  computations.  The  correlation  between  sec,  ship,  and  cable  stresses  is 
discussed  further  later  on  in  the  text  in  the  application  of  the  results  obtained  herein 
to  two  hypothetical  prototype  cases. 


RESULTS 


As  illustrated  above  in  the  theoretical  analysis,  the  parameters  influencing 
the  dynamic  stresses  can  be  tabulated  as  follows: 

Cable  Parameters 

L  =  maximum  length  of  cable 

S  =  material  cross-sectional  area  of  cable 
w  =  weight  of  cable  per  unit  length 
E  =  modulus  of  elasticity  of  cable 

=  allowable  maximum  dynamic  stress  in  the  cable 

Load  Parameters 

M  =  mass  of  the  load 
A  =  cross-sectional  area  of  the  load 

C  "  coefficient  of  mass 
m 

Cp  -  coefficient  of  drag 
p  =  density  of  sea  water 


Ship  or  floating  Platform  Motions  in  Heave 
|uo  -  amplitude  of  heove 

u>  =  frequency  of  heave  (rod/sec) 
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These  specific  variables  are  combined  as  follows: 
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Values  of  the  normalized  amplitude  of  the  maximum  dynamic  stress,  |  E^qxI, 
were  calculated  for  four  ranges  of  nondimensional  frequency,  u>\  The  ranges  were 
as  follows: 

tt/5  *  a:1  £  7.  Otr  in  increments  of  tt/5 
tt/10  s  w'  s  1 . 4tt  in  increments  of  tt/10 
0  s  k1  s  it/ 2  in  increments  of  tt/10 
0  s  uc1  s  0. 10  in  increments  of  0.01 


The  specific  values  of  M  used  in  the  computations  were  0.005,  0.01,  0.03, 
0.05,  0.10,  0.50,  1.00,  2.00,  and  5.  00.  For  low  ranges  of  a:',  additional  values 
of  M  were  used  as  shown  on  the  appropriate  graphs.  At  each  value  of  p,  calculations 
were  performed  for  fi  -  0.10,  0.30,  0.50,  0.70,  1.00,  3.00,  5.00,  and  7. 00  over 
the  two  higher  ranges  of  w',  and  values  of  0  =  0.25,  0.50,  1.00,  3.00,  5.00,  and 
7.00  over  the  two  lower  ranges  of  u;'.  The  results  of  these  computations  are  presented 
in  Figures  2  through  7  for  the  lowest  range  of  u;'  and  in  Figures  8  through  13  for  the 
second  lowest  range.  For  the  range  tt/10  *  u;'  s  1.4tt,  it  was  found  that  the  stress 
calculated  at  values  of  a;'  near  ft  were  zero.  Figure  14  illustrates  this  discontinuity. 
In  order  to  investigate  this  behavior,  Equations  7,  8,  9,  and  10  were  combined  and 
evaluated  as  u,'*  approached  nit,  where  n  -  1,  2,  3, . . .  Details  of  this  evaluation 
are  given  in  Appendix  B,  the  result  being  given  in  the  form 
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where  (9  is  the  parameter  defined  previously,  Equation  II,  and  K  =  nff/p.  The  »olid 
dot  in  Figure  14  is  the  stress  calculated  from  Equation  13  at  that  particular  value  of 
0  and  p.  The  values  of  |l‘  x|  as  determined  from  Equation  13  were  used  in  plotting 
Figures  15  through  22  as  ana  when  necessary.  The  apparent  discrepancies  in  the 
computer  calculations  are  considered  to  be  due  to  rounding-off  errors  inherent  in 
the  computational  procedures. 

Figure  23  is  included  to  indicate  the  variation  of  with  jS  for  a  particular 

value  of  M.  Figure  24  shows  the  variation  of  maximum  dynamic  stress  for  the  highest 
frequency  range  tr/5  s  a:'  <  7.  Off  for  particular  values  of  d  and  .  As  can  be  seen, 
the  maximum  dynamic  stress  is  highly  dependent  upon  the  nondimensional  frequency. 
Since  in  any  application  of  these  curves  the  peak  values  of  jIJnQX|  must  be  consid¬ 
ered,  the  results  for  this  range  of  frequencies  ore  presented  in  a  simplified  form  in 


Figures  25  through  32,  where  each  curve  is  drawn  through  the  maximum  values  cf 
l^maxl  'n  some  manner  as  the  dashed  line  in  Figure  24.  Each  graph  is  drawn  for 
various  values  of  U  at  a  given  0.  In  certain  instances,  it  was  again  found  necessary 
to  determine  at  frequencies  near  ff  and  its  multiples  by  use  of  Equation  13. 


DISCUSSION  OF  RESULTS 

The  results  obtained  from  the  computer  program  as  outlined  above  are  generally 
in  agreement  with  those  quoted  in  the  A.  D.  Little  report*  except  for  values  of  the 
normalized  maximum  dynamic  stress,  l^gxl'  corresponding  to  nondimensional 
frequencies,  u.*',  near  or  equal  to  3.  142.  The  computation  of  jl^xJ  a*  ~  n 
according  to  Equation  13  gave  results  which  compare  favorably  witn  those  obtained 
from  the  computer  program.  A  comparison  of  three  typical  results  is  given  in  Table  I. 

It  should  be  noted  that  the  value  of  j^gx)  computed  from  Equation  13  is  not 
necessarily  the  peak  value,  since  resonance  will  occur  at  nondimensional  frequencies 
other  than  ff  depending  upon  the  values  of  0  and  p.  When  p  approaches  zero,  the 
resonant  frequencies  approach  ff,  2ff,  3ff,  etc,  and  the  vclue  of  jE^xl  °*  determined 
from  Equation  13  may  then  be  interpreted  as  the  maximum  value.  This  con  be  seen  to 
be  true  by  inspection  of  Figures  15  through  22. 
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labia  I. 


Comparison  of  Typical  Values  of  jl^axj  °s  Computed  by  Digital 
Computer  Program  and  as  Calculated  From  Equation  13 


. .  . 

" 

P 

ii.'  I 

1  fT5ax' 

Computed  by  Digital 
Computer  Program 

|r  i 

.  1  ,maxi 
Cuicolaied  from 

Eqi*atior  13 

Percentage 

Difference 

0.50 

C.  10 

112  i  j  110.1 

1.80 

1.00 

0.  10 

142.32 

139.5 

2.01  i 

5,00 

0. 10 

514.09 

_ 

503.5 

2.  10 

— i 

The  parameters  jS  and  p  are  representative  of  the  damping  and  the  ratio  of  the 
weight  of  the  cable  to  the  virtual  mass  of  the  load  respectively.  As  was  noted  in  the 
A.  D„  Little  report,  the  variations  of  |l|noxj  with  a;',  /3,  and  U  are  in  agreement 
with  known  results  for  simpler  systems.  As  the  mass  of  the  load  is  decreased,  i.e., 
p  “*  ao,  the  system  reduces  to  that  of  a  free-ended  spring,  with  the  resonant  frequencies 
approaching  f/2  and  3tt/2,  etc.  As  th^  mass  is  increased,  p  ’•  0,  the  resonant  fre¬ 
quencies  approach  if,  2 IT,  etc.,  which  agrees  with  the  case  of  a  fixed-ended  spring. 

The  „amping  pnrameter,  0,  has  a  slight  effect  on  the  resonant  frequencies  but 
a  far  m  e  important  effect  or,  the  amplitude  of  j^cxj  at  resonance.  A  conclusion 
in  the  A.  D.  Little  report  indicated  that  the  maximum  dynamic  stress  amplitude  at 
res  lance  increases  when  the  damping  is  increased  beyond  a  certain  value.  From 
the  above  calculations  it  can  be  seen  that  the  amplitude  at  resonance  increases 
generally  with  increased  damping;  i.e. ,  there  is  no  minimum  amplitude  as  implied 
by  the  above  conclusion.  This  result  is  compatible  with  the  concept  that  as  the 
damping  increases,  the  system  becomes  equivalent  to  a  fixed-ended  spring  giving 
resonances  at  V,  2tt,  etc.,  and  amplitudes  tending  to  infinity,  restricted  only  by 
internal  and  external  damping  of  rhe  cable.  This  argument  considers  the  damping 
effects,  a  function  of  (3,  to  be  divorced  from  the  inertial  effects,  which  are 
dependent  u;*on  p. 

In  view  of  the  dependence  of  the  maximum  dynamic  stress  amplitude  on  ;S, 
and  since  0  depends  on  the  parameters  of  the  load  —  i.e.,  the  cross-sectional 
area,  the  mass,  and  the  density  of  sea  water,  which  are  fixed  —  and  on  the  coef¬ 
ficient  of  drag,  th,,  value  of  C[)  assumed  for  a  given  load  configuration  is  of 
particular  importance.  This  can  be  seen  from  the  results  given  above  where  a 
change  of  |3  from  1.00  to  3.00  results  in  a  change  in  |^ax|  at  resonance*  from 
142  to  370.  If  |^ax|  is  interpreted  as  an  allowable  stress  which  when  exceeded 
results  in  an  unsafe  condition  for  the  operation,  as  implied  in  the  design  procedure 
which  follows,  then  the  value  of  Cq  used  in  the  calculation  of  )3  becomes  critical. 

For  poor  hydrodynamic  shapes  such  as  blunt  bodies  or  open  frameworks,  it  is  not 
possible  within  the  present  state  of  development  of  theoretical  fluid  dynamics  to 
calculate  Cq  from  the  basic  equations  of  fluid  flow.  The  alternative,  therefore,  is 


to  resort  to  an  experimental  determination  of  the  coefficients  of  drag  for  the  particular 
load  configuration  in  question,  for  both  steady  and  oscillatory  motions.  There  exists 
little  experimental  data  on  the  appropriate  coefficients  of  drag  applicable  to  typical 
load  shapes  being  lowered  to  the  deep  ocean  floor.  A  series  of  experiments  directed 
toward  obtaining  such  data  thus  appears  justified. 

Similar  arguments  also  apply  to  the  values  assumed  for  the  coefficient  of  added 
mass,  Cm.  Summaries  of  information  pertinent  to  the  determination  of  the  drag  and 
added  mass  coefficients  are  presented  in  Appendixes  D  and  E  respectively. 

The  shape  of  each  curve  determine-  from  Equations  7,  8,  9,  and  10  differs 
from  these  obtained  for  linear  systems  in  that  a  second  peak  in  the  normalized  max¬ 
imum  dynamic  stress  occurs  at  nondimensional  frequencies  on  the  order  of  0. 10  to  0.40. 
The  significance  of  such  e  secondary  peak,  which  is  termed  herein  a  subharmonic 
response,  is  more  easily  discussed  in  relation  to  o  specific  design  example.  Two  such 
examples  are  given  below  following  a  proposed  design  procedure. 


PROPOSED  DESIGN  PROCEDURE 

The  parameters  required  in  the  design  procedure  are  those  applicable  to  the 
load  and  the  cable.  The'-  may  be  tabulated  as  follows: 

Load  Parameters 

M  =  mess  of  the  load 

A  =  cross-sectional  area  of  the  load  in  the  direction  of  motion  (ft^) 

Cp  =  coefficient  of  drag  applicable  to  the  load 

C  =  coefficient  of  added  mass  applicable  to  the  load 
m  r 

Cable  Parameters 

L  =  maximum  length  of  cable  (ft) 
max 

w  =  weight  of  cable  per  unit  length  (Ib/ft) 

E  =  modulus  of  elasticity  for  the  cable  (lb/in.^) 

I  ^  =  u^'m0*e  Ensile  strength  of  the  cable  (lb/in.^) 

F  -  safety  factor  for  maximum  operating  stress  in  the  cable 
p  =  density  of  sea  water  (Ib/ft^) 
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From  the  above  parameters,  the  following  may  be  determined: 


i 


(i) 


rd  ' 


ult  , 
p  ”  static 


where  equals  the  static  stress  in  the  cable,  including  that  due  to  the  cable 

itself  at  L  ,  and  equals  the  operational  maximum  dynamic  stress  allowable  with 
a  cable  of  ultimate  tensile  stress,  lyjj.,  and  a  safety  factor  of  F. 


the  velocity  of  sound  in  the  cable  in  feet  per  second;  and 


(3) 


4  CqP  A 
3tTc~M 

m 


a  constant.  Since  the  cable  length  varies  from  zero  to  and  because  the  design 
must  be  valid  for  all  lengths,  several  values  of  L  should  be  chosen  (L  =  L„  where 
n  =  1,  2,  3, . . . )  between  zero  and  Lmax.  Hence,  values  of  =  wL^C,,,  M  may 
be  calculated. 

Values  of  jU0|,  the  amplitude  of  the  cable  support-point  oscillation,  may  be 
selected  as  |UQ|  =  1.0,  2.0,  5.0,  10.0  feet.  A  table  may  then  be  set  up,  as 
illustrated  by  Table  II,  for  each  pn  corresponding  to  the  selected  Ln. 


Table  II.  Outline  of  Table  for  Computctior.  of  Relationship  Between 
Frequency  and  Amplitude  of  Oscillations,  Given  p,  L,  c, 
and  k 


{For  M  =  Mj,  L  =  l]t  c/L  =  c/Lj) 


Column  1 

Column  2 

Column  3 

Column  4 

Column  5 

m 

m 

wmmm 

(jl'1 

U)’  c 

“•'=  L 

1.0 

2.0 

5.0 

10.0 

15.0 
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The  use  of  the  table  is  as  follows.  Column  1  consists  of  the  amplitudes  of  the 
cable  support-point  oscillations  as  selected  above,  and  from  this.  Column  2  may  be 
calculated  by  multiplying  by  k  derived  previously.  Column  3  is  computed  for  the 
given  L,  Ij,  |UQ|,  and  E  values  appertaining  to  the  cable.  For  polypropylene  or 
nylon  cables,  it  may  not  be  possible  to  determine  and  E  directly.  Manufacturers' 
tables  generally  allow  the  computation  of  S  and  S  E,  where  S  is  the  material  cross 
section  of  the  cable.  Hence,  the  constant  ratio  Z<^/E  may  be  determined  and  used  in 
Column  3. 

Column  4  is  obtained  from  the  curves  of  Figures  2  to  7,  8  to  13,  15  to  22,  and 
25  to  32  for  the  particular  values  of  and  0  given  in  the  table.  Column  5  is  eval¬ 
uated  from  Column  4  for  the  particular  ratio  of  c/Ln .  Thus,  from  Columns  1  and  5  a 
curve  relating  the  allowable  amplitude  of  oscillation  to  the  frequency  of  that  oscilla¬ 
tion  may  be  drawn  for  the  particular  cable  length  used.  The  process  may  then  be 
repeated  for  other  cable  lengths.  A  judicious  choice  of  cable  length  can  serve  to 
reduce  the  numerical  calculation  to  a  minimum. 

It  may  be  assumed  that  operating  conditions  lying  on  or  below  the  curve  are 
safe,  with  the  safety  factor,  F,  as  defined,  and  that  operating  conditions  lying  above, 
to  the  right  of,  the  curve  are  unsafe. 

In  line  with  comments  raised  in  the  Discussion  of  Results,  the  unknown  parameters 
are  Cp  and  Cm.  With  existing  deficiencies  •  «  data  giving  Cp  and  Cm  for  various  load 
configurations,  they  must  of  necessity  be  estimated.  See  Appendixes  D  and  E. 

If  it  should  occur  in  a  design  problem  that  the  ranges  of  u:1,  0,  or  p  are  not 
covered  in  the  graphs  developed  in  this  study,  then  the  appropriate  curves  may  be 
calculated  for  specific  values,  or  ranges  of  values,  of  those  parameters  by  use  of 
the  digital  computer  program  used  in  obtaining  the  results  quoted  herein.  Details 
of  this  program  are  given  in  Appendix  C. 

Two  examples  of  the  application  of  the  above  design  procedure  are  given  below. 


APPLICATION  OF  PROPOSED  DESIGN  PROCEDURE  TO  TWO  PROTOTYPE 
EXAMPLES 

The  application  of  the  proposed  design  procedure  to  two  hypothetical  prototype 
examples  is  demonstrated  below  for  polypropylene  and  steel  cable  respectively.  The 
parameters  used  in  these  examples  are  such  as  to  enable  the  design  procedure  to  be 
carried  out  using  the  curves  presented  previously  in  the  results. 

Design  Example  Using  Polypropylene  Cable 

The  parameters  of  the  load  are  given  as 
M  =  5.0  tons  =  10,000  lb 
A  =  12  ft  x  12  ft  =  144  ft2 
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The  cable  parameters  chosen  are 


L 

max 

w 

SE 

szd 

c 


k 


20, 000  ft 
0. 90  Ib/ft 
240, 000  lb 
10,000  lb 


4  CDP  A 

3"CmM 


0.50 


2, 930  ft/sec 


and,  therefore,  jS=0.50|Uo|. 

The  cable  data  used  herein  was  obtained  from  the  August  1964  "Braided  Rope 
and  Cordage  Catalog"  of  the  Samson  Cordage  Works,  Boston,  Mass.  It  is  now 
convenient  to  select  cable  lengths  such  that  values  of  where  ji  =  wL/Ma , 
coincide  with  those  used  in  deriving  the  curves  presented  previously  in  the  results. 
These  values  are  given  in  Table  III. 


Table  III.  Cable  Lengths  Used  in  Design  Example  for 
Polypropylene  Cable 


L  (ft) 

J£ 

cA 

16, 660 

1.00 

0.1758 

8, 330 

0.50 

0.3516 

1,660 

0. 10 

1.7580 

833 

0.05 

3.5160 

166 

0.01 

17.5800 

It  is  now  possible  to  set  up  Table  IV  corresponding  to  Table  II  given  in  the 
proposed  design  procedure  by  specifying  that  input  amplitudes  of  oscillation,  |UQj, 
of  1.0,  2.0,  6.0,  10.0,  and  14.0  feet  will  be  considered.  Knowing  that  3  =  k  |U0|, 
Column  2  of  Table  IV  may  be  calculated.  The  volues  of  | Z^!nax  |  =  L  |U0|  E  can 
be  derived  for  appropriate  values  of  l  and  |U0|  as  shown  in  Column  3  of  the  table. 
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* 


From  inspection  of  the  table,  it  can  be  seen  that  in  view  of  the  relationship  between 
the  various  cable  lengths  chosen,  values  of  jl^ax|  at  these  different  lengths  are  quite 
simply  related  —  thus  facilitating  the  design  calculations. 

Values  of  the  nondimensional  frequency,  u;',  are  then  entered  in  Column  4 
by  the  use  of  Figures  2  to  13,  15  to  22,  and  25  to  32.  That  is,  for  a  particular  p 
and  d/  the  value  of  a.'1  corresponding  to  ma/  ^  found.  Hence,  the  circular 

frequency,  a.',  can  be  calculated  and  entered  in  Column  5.  From  the  completed 
table,  the  relationship  as  a  function  of  cable  length,  can  be  drawn  between  input 
amplitude  of  oscillation,  )  U  ] ,  and  the  allowable  circular  frequency,  uc,  of  that 
amplitude  —  i.e.,  the  circular  frequency  at  which  the  oscillation  can  occur  such 
the  the  maximum  dynamic  stress  in  the  cable  is  less  than  or  equal  to  the  design 
dynamic  stress.  Figure  33  shows  this  relationship  for  the  computation  given  in 
Table  IV.  The  significance  of  these  results  is  discussed  below  together  with  that  of 
the  following  design  example. 

Design  Example  Using  Steel  Cable 


The  relevant  load  parameters  are  given  as  follows: 
M  =  20  tons  =  40, 000  lb 
A  =  600  ft2 
CD  =  2.0 
C  =  1-5 

m 

The  appropriate  cable  parameters  were  chosen  to  be 
Lmax  =  20,000  ft 
pc  =  550  Ib/ft3 
W  =  7. 64  Ib/ft 
E  =  15  x  106  psi 
Id  -  40,000  psi 
c  --  11,200  ft/sec 


k 


4CDpA 
3iTCm  M 


0.50 


and,  therefore,  &  ~  0.50  j'J^ 
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Table  IV.  Design  Calculations  for  10,000-Pound  Load  to  20,000  Feet 
Using  Polypropylene  Cable 
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7.00  2.48  0.68tr  7.511  0.496  0.0280  0.492  0.248  0.0130  0.457 


As  in  the  first  example,  cable  lengths  are  chosen  to  give  values  of  n  coincident 
with  those  used  in  the  calculation  of  Figures  2  to  13,  15  to  22,  and  25  to  32.  These 
lengths  and  the  corresponding  values  of  H  and  c/L  are  given  in  Table  V. 


Table  V. 

Cable  Lengths  Used  in  Design  Example 
for  Steel  Cable 

L  (ft) 

c/L 

15,700.0 

2.000 

0.713 

7, 850.0 

1.000 

1.426 

3, 925.0 

0.500 

2.850 

1,963.0 

0. 250 

5.600 

785.0 

0. 100 

14.260 

393.0 

0.050 

28.500 

236.0 

0.030 

47.500 

78.5 

0.010 

142.600 

39.3 

0. 005 

285.000 

A  similar  table  to  that  derived  in  the  previous  design  example  may  now  be 
set  up  and  the  numerical  computations  performed  as  above.  These  calculations  are 
summarized  in  Table  IV,  and  the  relationship  between  the  input  amplitude  of 
oscillation,  [U0|,  and  the  allowable  circular  frequency,  to,  of  that  amplitude  is 
illustrated  in  Figure  34. 


DISCUSSION  OF  RESULTS  OF  THE  APPLICATION  OF  THE  PROPOSED 
DESIGN  PROCEDURE 

Figures  33  and  34  show  the  results  obtained  in  the  application  of  the  proposed 
design  procedure  to  the  two  hypothetical  cases  described  above. 

In  any  load-lowering  operation,  the  cable  used  will  have  a  certain  known 
ultimate  load  at  which  the  cable  could  be  expected  to  break.  With  repeated  use 
of  a  particular  cable,  this  ultimate  load  will  decrease.  Hence,  for  any  lowering 
operation  a  safety  factor  must  be  chosen  defining  a  load,  or  corresponding  stress, 
which  should  not  be  exceeded.  This  allowable  working  stress  is  assumed  to  include 
static  stress  due  to  both  the  load  and  cable  —  the  latter  being  negative  in  the  case 
of  a  buoyant  cable  such  as  polypropylene  —  arid  the  dynamic  stress. 

This  discussion  and  the  design  procedure  proposed  earlier  in  the  report  are 
based  on  the  assumption  that  a  particular  maximum  dynamic  stress  is  given  which 
should  not  be  exceeded  during  the  lowering  or  raising  operation. 
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Figures  33  and  34  indicate  the  allowable  frequency  of  an  input  oscillation  at 
various  amplitudes  of  this  oscillation  as  a  function  of  the  cable  length,  such  that  the 
stipulated  maximum  dynam'c  stress  is  not  exceeded.  Alternatively,  they  indicate 
the  maximum  amplitude  of  an  input  oscillation  at  a  particular  frequency  such  that  a 
maximum  dynamic  stress  is  not  exceeded. 

One  of  the  most  important  problems  in  the  interpretation  of  these  graphs  is  the 
selection  of  appropriate  input  conditions,  |U0|  and  to,  corresponding  to  the  response 
of  the  vessel  used  for  the  lowering  or  raising  operation  of  the  sea  state  in  which  the 
operation  is  carried  out. 

It  is  apparent  that  the  frequency  range  given  in  Figures  33  and  34  is  far  greater 
than  that  which  would  be  expected  under  operational  conditions.  Data  is  available 
in  the  literature  (Kaplan  and  Putz,  2  Pierson  and  Holmes^)  on  the  response  of  the 
Cuss*-!  drilling  barge  to  various  sea  states.  In  sea  state  5,  the  range  of  frequencies 
of  oscillation  in  heave  is  0.40  to  1.40  radians  per  second  with  root -mean -square 
values  of  heave  of  1. 8  and  1.3  feet  at  headings  of  90  degrees  and  0  degrees  to  the 
wind  respectively.  These  values  imply  expected  maximum  amplitudes  of  oscillation 
of  7.2  and  5.2  feet  during  a  4-hour  period  on  station  (Pierson  and  Holmes).  It  should 
be  noted,  however,  that  combinations  of  heave  and  roll  oscillations  could  easily  pro¬ 
duce  oscillations  greater  than  this  if  the  load-lowering  operation  is  carried  out  using 
the  boom  over  the  s»de  of  the  vessel.  The  above  values  are  used  here  to  Illustrate 
the  interpretation  of  Figures  33  and  34. 

Referring  to  Figure  33  for  the  design  example  using  a  polypropylene  cable, 
it  can  be  seen  that  for  a  heave  amplitude  of  7.  2  feet  at  frequencies  of  1.40  and 
0.40  radians  per  second  —  to  cover  the  entire  frequency  range  —  the  design  dynamic 
stress  will  be  exceeded  at  a  cable  length  of  less  than  200  feet  for  cc  =  1.40.  For 
frequencies  less  than  1.40  radians  per  second  down  to  0.40  radian  per  second,  it  is 
estimated  from  the  curves  that  the  design  dynamic  stress  will  be  exceeded  at  cable 
lengths  which  gradually  decrease  from  200  feet.  Thus  the  operation  will  be  unsafe 
relative  to  the  prescribed  maximum  dynamic  stress  for  a  cable  length  of  less  than 
200  feet.  This  condition  has  been  fully  recognized  in  the  design  of  various  lowering 
operations  conducted  by  this  laboratory.  The  dynamic  stress  will  then  be  less  than 
the  design  stress  until  the  length  of  the  cable  reaches  a  value  of  approximately 
5,000  feet.  According  to  the  results  obtained  above  the  operation  will  then  become 
unsafe,  if  the  frequency  is  1.40  radians  per  second  and  the  amplitude  is  7.  2  feet.  * 
For  a  frequency  of  1.00  radian  per  second  the  operation  becomes  unsafe  at  a  depth 
of  7,000  feet.  The  reasons  for  this  result  are  difficult  to  visuclize,  but  may  be 
explained  from  both  the  mathematical  and  physical  points  of  view.  In  terms  of  the 
numerical  computations  carried  out  as  shown  in  Table  VI,  the  regression  of  the  curves 
given  in  Figures  33  and  34  is  due  to  the  fact  thot  when  the  maximum  dynamic  stress, 

*  These  comments  ore  based  on  an  approximate  interpolation  beiween  the  curve 
labeled  6  feet  and  10  feet  in  Figure  33  for  |  U0|  ~  7.2  feet. 
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I VI'  increases,  say  doubles,  due  to  an  increase  in  the  cable  length  of  L  to  twice 
L,  the  appropriate  nondimensional  frequency,  a:1,  at  the  corresponding  value  of  n 
is  not  twice  the  a,'1  at  the  shorter  length.  This  results  in  a  lower  circular  frequency, 
a;,  at  the  longer  cable  length.  This  applies  to  the  larger  values  of  {correspond¬ 

ing  to  these  lengths.  A  similar  argument  applies  at  cable  lengths  lower  than  200  feet. 
For  intermediate  cable  lengths  a  relatively  small  variation  in  |  T^qxJ  results  in  a 
significant  increase  in  cc1,  resulting  in  an  increase  in  the  circular  frequency,  u:. 

This  applies  to  a  range  of  a.'1  from  0.4tr  to  0. 8ff  —  i.e.,  to  values  of  |Z^ax|  larger 
than  thore  corresponding  to  the  first  peak  in  Figures  2  to  13  and  15  to  22,  but  less 
than  those  associated  with  the  peak  responses  occurring  at  a;1  approximately  equal 
to  3. 142.  As  oc1  values  tend  toward  inis  value,  significant  changes  in  j  T^qx]  result 
in  larger  variations  in  circular  frequency  for  different  cable  lengths  due  to  the  influ¬ 
ence  of  the  c/L  ratio  in  computing  the  latter.  Hence  the  variation  of  cable  stress 
with  cable  length  depends  upon  the  shapes  of  the  computed  curves  relating  j Z*  I 
tow  . 

The  physical  behavior  of  the  cable  assembly  may  be  described  by  reviewing 
the  significant  results  obtained  by  Little  1  and  by  Whicker. 4  The  latter  is  a  rather 
simplified  (i.e.,  no  damping)  theoretical  analysis  of  the  effect  of  ship  motion  on 
mooring  cables  in  deep  water.  On  Figure  33,  values  of  ui'  corresponding  to  n/2  and 
V  are  indicated,  as  well  as  the  roots  of  the  equations  tantu1  =  \i/x'  and  tanoc'  = 

The  case  of  the  fixed-ended  spring  for  zero  damping  (J  =  0. 0)  is  shown  by  Whicker 
and  by  Little  to  have  resonant  frequencies  of  a;'  =  n,  2 it, .  ,.nv  when  the  relative 
mass,  fi,  of  the  cable  and  payload  is  decreased  infinitely  (he**  the  payload  mass 
increases  indefinitely).  The  case  of  the  free-ended  spring  for  zero  damping  (8-  0.0) 
is  represented  by: 

1.  Values  of  a:'  equal  to  *r/2,  3it/2. .  .  ntt/2  when  the  relative  mass  is  increased 
indefinitely.  When  the  damping  8  =  0. 0,  it  is  shown  by  Little  that  values  of 
tx'  equal  to  it/ 2,  3ff/2.  ..nTt/2  result  in  infinite  dynamic  stress  for  infinitely 
large  values  of 

2.  Values  of  a.'1  corresponding  to  the  roots  of  the  equation  tan  u,''  =  -n/x'. 
Whicker  shows  that  for  finite  values  of  4,  the  least  root  of  the  above  equa¬ 
tion  results  in  infinite  total  stress. 

In  oddition.  Little  shows  that  for  finite  values  of  n,  and  for  zero  damping, 
infinite  total  stress  will  result  when  the  roots  of  the  equation  tan  u,'1  =  y/<x'  ore 
satisfied.  This  is  in  contrast  to  the  results  given  by  Whicker,  yet  both  formulations 
appear  to  be  correct.  Thus,  it  would  appear  that  the  Little  ond  Whicker  results  are 
compatible  for  long  coble  lengths,  but  ttat  the  Whicker  results  are  not  applicable 
for  short  lengths,  since  all  of  recorded  data  supports  Little's  conclusions.  In  any 
event,  the  least  of  the  roots  of  the  above  equation  lies  between  0  and  ff/2.  As  ji  is 
increased  indefinitely,  the  roots  of  the  above  equation  approach  it/ 2,  3ir/2. . .  nft/2. 
There  is  no  comparable  analogy  with  the  well-known  results  of  simpler  systems. 
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Table  VI.  Design  Calculations  for  40,000-Pound  Load  to  20,000  Feet 
Using  Steel  Cable 


0. 50  0.6276  0.1020  4.  WO  0.2093  0.0385  5.490  0.1046  0.0205  5.840 

1.00  0.3138  0.0715  3.396  0.1046  0.0255  3.636  0.0523  0.0125  3.560 

3.00  0.1044  0.0305  1.448  0.0349  0.0100  1.426  0.0174  0.0030  0.854 

5.00  0.0628  0.0190  0.903  0.0209  0.0050  0.713  0.0105  0.0019  0.542 

7.00  0.0448  0.0140  0.665  0.0149  0.0040  0.570  0.0075  0.0015  0.428 


Beginning  with  large  lengths,  the  cable  assembly  behaves  as  a  fixed-ended 
spring,  and  the  role  of  damping,  which  is  dependent  on  the  input  amplitudes,  is  to 
provide  an  additional  margin  of  safety.  Thus  higher  frequencies  for  the  same  length 
at  constant  input  amplitudes  can  be  tolerated  without  danger  of  breaking  the  cable. 

As  the  length  decreases,  a  transition  occurs  and  the  cable  assembly  takes  on 
the  characteristics  of  a  free-ended  spring.  Again  the  role  of  damping  is  to  provide 
an  additional  margin  of  safety  against  failure. 

Finally  at  very  short  cable  lengths,  the  cable  appears  to  behave  liko  a  rigid 
connection  between  source  and  load  rather  than  as  a  "spring. "  In  this  cas<»,  the 
role  of  damping  (drag)  is  reversed  in  that,  for  constant  input  amplitudes,  smaller 
frequencies  car  be  tolerated  than  for  the  case  when  j3  =  0. 0. 

All  of  the  above  trends  are  confirmed  by  the  calculations  for  both  the  damped 
and  undamped  cases  as  the  length  L  approaches  zero,  as  indicated  by  Figure  33. 

For  comparative  purposes,  a  few  results  obtained  by  applying  the  method  developed 
by  Whicker  for  a  free-ended  cable  are  included.  The  mathematical  formulation  by 
Whicker  does  not  include  the  damping  term.  Thus  a  direct  comparison  of  the  effect 
of  damping  is  available. 

It  is  recognized  that  the  analysis  used  for  this  investigation  —  that  given  in 
the  report  by  A.  D.  Little,  Inc.  —  does  not  accurately  describe  the  prototype  situ¬ 
ation  by  virtue  of  the  linearization  reauired  in  the  drag  term  j  :'u/ot'|  (c'u/^t'),  which 
is  necessary  in  order  to  solve  the  basic  equation  of  motion.  In  the  present  state  of 
the  art  concerning  nonlinearly  damped  oscillations  there  appear,  to  be  no  alternative 
to  the  linearization.  Further  theoretical  analysis  is  considered  necessary  and  justified 
in  order  to  resolve  the  question  of  safety  of  lowering  systems  as  ihe  length  of  the 
cable  increases.  Such  an  analysis  would  consider  the  use  of  analogue  computations 
to  allow  retention  of  the  nonlinear  |  ?u/  't*  |  (cu/  '*')  term. 

Similar  interpretations  regarding  the  safety  of  a  load-lowering  operation  may 
well  apply  at  greater  depths  for  lower  frequencies  if  the  curves  presented  in  Figure  33 
were  extended.  From  the  general  shape  of  the  curve  shown  and  for  a  given  input 
amplitude  oscillation  there  is  a  limiting  input  frequency  below  which  the  operation 
is  safe  (the  maximum  dynamic  stress  is  not  exceeded)  down  to  a  certain  depth.  For 
a  particular  lowering  operation  it  may  be  possible  to  use  a  working  vessel  which  has 
little  or  no  response  to  excitations  above  this  limiting  frequency.  Alternatively,  the 
operation  may  be  carried  out  in  a  lower  sea,  but  this  does  not  necessarily  imply  a 
maximum  frequency  of  input  oscillation,  althotgh  it  does  imply  a  diminished  ampli¬ 
tude  of  oscillation  which  would  render  the  operation  safe. 

A  similar  discussion  may  also  be  applied  to  Figure  34  for  the  cose  of  a  steel 
coble,  although  an  unsafe  condition  is  not  likely  to  occur  until  depths  of  20,000  feet 
even  with  input  amplitudes  of  14  feet.  In  this  case,  however,  an  additional  factor 
must  be  considered,  nomely  the  very  significant  increase  in  static  stress  due  to  the 
weight  of  the  cable.  As  far  as  the  dynamic  stress  is  concerned,  it  moy  be  concluded 
that  the  operation  moy  be  safe  for  coble  lengths  greater  than  100  feet  when  input 
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amplitudes  are  less  than  6  feet  at  a  maximum  frequency  of  1.40  radians  per  second. 

As  the  cable  length  diminishes  to  zero  it  would  appear  that  the  design  dynamic  stress 
will  be  exceeded  by  an  increasingly  large  amount.  From  Figures  33  and  34  it  is  seen 
that  for  both  examples  at  shorter  lengths  of  cable,  a  greater  amplitude  of  oscillation 
is  relatively  less  safe  than  a  small  amplitude,  as  expected.  Further  calculations  are 
required  to  determine  the  form  of  the  curves  as  the  cable  length  approaches  zero. 

On  13  April  1965,  a  Submersible  Test  Unit  (STU)  ioaded  with  racks  of 
specimens  was  lowered  to  a  depth  cf  2,500  feet  by  the  Deep  Ocean  Engineering 
Division  of  NCEL.  The  record  of  cable  tensions  during  the  lowering  operation  h 
shown  on  Figure  35.  The  weight  of  the  load  in  water  wes  approximately  5,500  pounds, 
the  structure  being  in  the  form  of  an  open  truss,  its  base  consisting  of  two  fiat  plates 
with  a  total  cross-sectional  area  of  approximately  150  square  feet.  The  load  Vvjs 
lowered  on  a  1 . 3-inch-diameter  polypropylene  cable.  From  a  depth  of  450  feet  the 
descent  was  carried  out  at  a  steady  rate  of  132  feet  per  minute.  At  depths  shallower 
than  450  feet  the  lowering  operation  was  intermittent.  The  wave  excitation  was 
estimated  to  be  in  the  form  of  a  6-  to  10-foot  swell  with  periods  of  10  to  12  seconds. 

A  brief  discussion  of  various  properties  of  the  record  in  the  light  of  the  calculations 
given  above  is  pertinent.  The  parameters  of  the  cable  and  load  were  approximately 
those  used  for  the  first  design  example,  although  values  of  the  coefficients  of  drag 
and  virtual  mass  do  not  necessarily  agree. 

The  record  shows  an  expected  decrease  of  dynamic  stress  with  increasing  depth. 
The  immediate  reduction  in  the  mean  load  at  point  B  of  approximately  2,000  pounds 
corresponds  to  a  drag  force  on  the  structure  with  a  coefficient  of  drag  equal  to  2.76 
at  a  vertical  velocity  of  132  feet  per  minute.  Since  the  load  consisted  of  racks  of 
specimens  orientated  pa-a!lel  to  the  flow,  this  drag  coefficient  is  not  necessarily  that 
corresponding  to  form  drag  alone;  however,  the  contribution  of  tangential  drag  on 
the  cable  can  be  shown  to  be  negligible.  The  steady  reduction  of  mean  *oad  from 
450  to  2,500  (eet  is  due  to  the  gradual  removal  of  the  weight  of  a  2,500-foot-long 
steel  cable  suspended  from  the  base  of  the  STU  to  the  ocean  floor. 

From  the  design  cxampie  given  previously  the  expected  maximum  dynamic 
stresses  in  the  lowering  cable  were  calculated  at  cable  lengths  of  83,  830,  and 
1,660  feet  assuming  that  |UQ  j  equals  6.0  feet.  Table  VII  summarizes  these  compu¬ 
tations  which  were  carried  out  for  each  cable  length  by  determining  the  period  of 
the  cyclic  stress  from  the  record.  Values  cf  jl^0x|  were  then  found  from  the  curves 
given  In  the  results,  for  ihe  appropriate  d  and  values.  Hence,  the  maximum 
dynamic  stresses  were  determined  as  shown  in  Table  VII  and  superposed  on  the  stress 
record  os  given  in  Figure  35. 

In  view  of  the  uncertainty  with  respect  to  values  of  Cp,  Cm,  and  |U0|  for  the 
operation  represented  by  Figure  35,  the  comparison  of  the  calculated  dynamic  stresses 
with  those  determined  experimentally  is  considered  to  be  fair.  It  is  noted  that  at 
830  and  1,660  feet  the  stresses  recorded  exceed  those  calculated  by  a  factor  of  2. 

As  the  load  is  lowered,  only  one  variable  changes,  that  being  the  cable  length. 
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However,  the  frequency  and  amplitude  of  the  cyclic  stress  change  with  cable  length. 
This  may  well  result  in  a  significant  change  in  the  coefficient  of  drag  and  mass  since 
the  amplitude  and  velocities  of  vertical  oscillation  will  vary.  This  is  possibly  one 
cause  of  the  discrepancies  noted. 


Table  VII.  Summary  at  Calculations  for  STU  Lowering  Operation 
to  2,500  Feet 


Cable  Lengfh,  L 

83.  0  ft 

830. 0  ft 

1,660.0  ft 

Parameter, 

0.005 

0.05 

0. 10 

Parameter,  j8 

3.00 

3.00 

3.00 

c/L 

0.  3283 

0.283 

0.566 

Period  of  Oscillation,  T 

5.45  sec 

8. 00  sec 

9. 23  sec 

Frequency  of  Oscillation,  oo 

1.  153  rads/sec 

0. 786  rad/sec 

0.681  rad/sec 

i  a)  L 

a:  -  — 
c 

0.0326 

0.222 

0.385 

I'  From  Curves 

1  max 

0.660 

1.165 

1.225 

S  -  Dynamic  Load 

±3,752.0  lb 

±707.4  lb 

±371.9  lb 

Static  Load 

5,500  lb 

3,400  lb 

3, 100  lb 

Static  +  Max  Dynamic  Load 

9,252  lb 

4,  107  lb 

3,472  lb 

Static  -  Max  Dynamic  Load 

1,748  lb 

2,693  lb 

2, 728  lb 

From  the  record  given  in  Figure  35  there  also  appears  to  be  a  somewhat 
irregular,  long-pericd  reinforcement  of  the  dynamic  stress  amplitude.  During  the 
lowering  operation  it  was  not  possible  to  record  the  actual  motions  of  the  load,  and 
thus  any  departure  from  purely  vertical  motions  are  unknown.  At  a  lowering  velocity 
of  132  feet  per  minute  it  is  quite  possible  that  the  load  motion  will  become  unstable 
and  that  sidewise  oscillations  will  occur,  giving  rise  to  complex  lift  forces  on  the 
structure  with  resulting  variations  in  the  cable  tensions.  There  is  no  evidence  in 
Figure  35  to  support  the  result  shown  in  Figure  34  that  dynamic  stresses  may  increase 
at  greater  depths,  because  of  the  limited  depth  of  2,500  feet  to  which  the  STU  was 
lowered.  There  is  a  pressing  need  for  a  load-lowering  operation  to  be  carried  out  in 
which  both  the  cable  tensions  and  the  load  motions  are  recorded,  thus  allowing  a 
more  definite  analysis  of  the  results.  Preferably  this  full-scale  experiment  would  be 
carried  out  with  a  formalized  body  shape  so  that  coefficients  of  drag  and  mass  could 
be  more  accurately  estimated. 
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FINDINGS 


As  a  result  of  the  investigation  described  above,  it  was  found  that: 

1.  The  theoretical  analysis  presented  by  Arthur  D.  Little,  Inc.,  in  Project  Trident 
Technical  Report  No.  1370863  of  the  Bureau  of  Ships1  was  readily  adaptable  to  a 
design  procedure  for  heavily  loaded  cables  for  deep  ocean  emplacement  or  recovery 
operations. 

2.  This  design  procedure  is  relatively  straightforward,  but  the  results  of  its  applica¬ 
tion  require  confirmation  by  prototype  measurements. 

3.  The  most  important  unknowns  in  the  input  parameters  for  the  design  procedure  are 
the  coefficients  of  drag  and  virtual  mass,  which  at  this  stage  must  of  necessity  be 
estimates,  again  requiring  confirmation  by  prototype  measurements. 

4.  For  the  two  prototype  examples  given,  the  operations  are  unsafe  with  respect  to 
the  specified  maximum  allowable  dynamic  stresses  when  the  load  is  between  zero  and 
approximately  200  feet  below  the  ocean  surface. 

5.  Using  polypropylene  cable,  there  exists  a  possibility  that  under  certain  conditions 
of  motion  of  the  working  vessel  the  operation  will  again  become  unsafe  at  depths  on 
the  order  of  5,000  feet  and  greater.  The  reason  is  not  difficult  to  visualize.  As  the 
length  increases  the  natural  period  also  increases,  eventually  corresponding  to  periods 
of  the  exciting  waves.  If  the  damping  is  small,  very  large  stresses  may  be  induced  by 
relatively  small  input  amplitudes. 


CONCLUSIONS 

It  is  concluded  that  the  proposed  design  procedure  is  applicable  to  heavily 
loaded  cables  in  deep  ocean  emplacement  or  recovery  operations.  However,  the 
results  of  its  application  may  not  be  considered  rigorous  in  view  of  the  estimates  of 
values  of  coefficients  of  drag  and  mass  required  for  the  calculations.  In  order  to 
make  the  design  procedure  applicable  to  a  prototype  situation  with  a  greater  degree 
of  confidence,  it  is  considered  necessary  to  make  measurements  of  coble  tensions 
and  load  and  ship  motions  during  a  full-scale  operation  to  provide  a  basis  for  com¬ 
parison  between  theory  and  prototype. 
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Figure  1.  Definition  sketch  cf  lowering  or  raising  operation. 
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Figure  2.  Variation  of  normalized  maximum  dynamic  stress,  |^ax|/  with 
nondimensional  frequency,  a)1,  for  /3  =  0.25  over  the  range 
0  £  a)'  *  0.10  for  values  of  v  as  indicated. 

24 


i 


•i' 

t 


Jf 


i 


Normalized  Maximum  Dynamic 


lit 


♦■tj 


♦Mi;' 


m-t *■) 


Figure  6.  Voriotioo  of  normalized  maximum  dynamic  stress,  |^!IKJX|,  with 
nondimereional  frequency,  u>',  for  0  =  5.00  over  the  range 
0  *  UJ*  *  0.10  for  values  of  u  as  indicated. 


Figure  7.  Variation  of  normalized  maximum  dynamic  stmts,  I^qxJ/  with 
nondimemional  frequency,  a;',  for  fi  =  7.00  over  the  range 
0  W  s  0.10  for  values  of  M  o*  indicated. 
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Normalized  Maximum  Dynamic  Stress,  -  ^ , 


NondimenttoneJ  Frequency.  *  \  £  • 

Figure  8.  Variation  of  normalized  maximum  dynamic  stress,  i^maxl' 
nondimensional  frequency,  for  8  =  0.25  over  the  range 
O4*'1  s.  t/2  for  values  of  ji  as  indicated. 
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Figure  10.  Variation  of  normalized  maximum  dynamic  stress,  IlL  I.  with 

_ i? _ ? _ i  r  i  .  a  .  ..  i  rnaxi' 
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0sO!'  s  it/ 2  for  values  of  as  indicated. 
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Normalized  Maximum  Dynamic  Stress 


M  =  0.04 
u  =  0.03 

M  =  0.02 
M  -  0.01 
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Nondimensional  Frequency,  ■  (— ~-) 

Figure  11.  Variation  of  normalized  maximum  dynamic  stress,  |^axj/  with 
nondimensional  frequency,  ic',  for  /3=  3.00  over  the  range 
0  s  co*  *  ff/2  for  values  of  ji  as  indicated. 
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Figure  13.  Variation  of  normalized  maximum  dynamic  stress,  I^qxI#  with 
nondimensiona!  frequency,  a;',  for  fi  =  7.00  over  the  range 
0  *  w'  s  ff/2  for  values  of  ja  as  indicated. 
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Figure  15. 
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Variation  of  normalized  maximum  dynamic  stress.  1 1*  . 

•  max  w 

nondimensional  frequency,  a;',  for  0  -  0.10  over  the  range 
tt/10  *  O)1  s  1.4tr  for  values  of  4  as  indicated. 
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Figure  17.  Variation  of  normalized  maximum  dynamic  stress,  j^nax|^ 
nondimensional  frequency,  u for  0  =  0.50  over  the  range 
tr/10  *  a>'  *  1.4tr  for  values  of  u  as  indicated. 
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Figure  18.  Variation  of  normalized  maximum  dynamic  stress,  ,  with 

nondimensional  frequency,  u;1,  for  fl  =  0.70  over  the  range 
tr/10  s  a.'1  *  1.4tr  for  values  of  n  as  indicated. 
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Figure  20.  Voriation  of  normalized  maximum  dynamic  stress,  , 

nondimensional  frequency,  u;',  for  $  ~  3.00  over  the  range 
tr/10  *■  a;'  i  1.4#  for  volues  of  u  os  indicated. 
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Figure  25.  Variation  of  normalized  maximum  dynamic  stress,  with 

nondimensional  frequency,  w1,  for  £  -  0.10  over  the  range 
it/ 5  *  u >'  *  7.0ff  for  values  of  4  at  Indicated. 
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Figure  26. 
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Variation  of  normalized  maximum  dynamic  stress,  j^axl'  with 
nond intentional  frequency,  w1,  for  /3  =  0.30  over  the  range 
if/ 5  *  w*  *  7.0ff  for  valuet  of  H  as  indicated. 
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Nondim«ntionol  Fr«qu«ncy,  (*>  *  \ 

Figure  27.  Variation  of  normalized  maximum  dynamic  stress,  |^ax  , 
nondimensional  frequency,  go',  for  0  =  0.50  over  the  range 
ff/5  *  co'  *  7.0ff  for  values  of  ji  as  indicated. 
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Figure  28,  Variation  of  normalized  maximum  dynamic  stress, 

nondimensional  frequency,  a;',  for  0  ~  0.70  over  tne  i 
tt/5  4  a)*  4  7. Off  for  values  of  u  as  'ndicated. 
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Figure  29.  Variation  of  normalized  maximum  dynamic  stress,  |^naxj , 
nondimensiona I  frequency,  a.'1,  for  0  =  l.OGove  the  range 
ff/5  s  uJ*  *  7.0f r  for  values  of  4  as  indicated. 
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Figure  30.  Variation  of  normalized  maximum  dynamic  stress,  U^x)  #  with 
nondimensional  frequency,  ui',  for  $  =  3.00  over  the  range 
tr/5  ^  w*  s  7. On  for  values  of  n  as  indicated. 
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Figure  31.  Variation  of  normalized  maximum  dynamic  stress,  |l^naxj ,  with 
nondimensional  frequency,  u;',  for  0  =  5.00  over  the  range 
ff/5  4  0;'  4  7. Of  for  values  of  4  as  indicated. 
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Figure  32.  Variation  of  normalized  mcximum  dynamic  stress, 
nondimensional  frequency,  uc',  for  0  =  7.00  over 
ff/5  *  a?’  *  7. Off  for  values  of  u  as  indicated. 
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Figure  34.  Results  for  design  example  using  steel  cable 


Figure  35.  Record  of  cable  tension  during  5TU  lowering  operation. 

Comparison  with  calculated  maximum  dynamic  stresses. 


vring  STU  lowering  operation, 
•d  maximum  dynamic  stresses. 


Appendix  A 


ANALYSIS  OF  CABLE  AND  PAYLOAD  DYNAMICS  FOR 
RAISING  AND  LOWERING  LOADS  IN  THE  DEEP  OCEAN 


The  analysis  presented  here  is  essentially  that  given  in  BuShips  Technical 
Report  No.  1370863,  "Stress  Analysis  of  Ship-Suspended  Heavily  Loaded  Cables  for 
Deep  Underwater  Emplacements,"  by  Arthur  D.  Little,  Inc.  ^ 

The  cable  and  load  are  considered  to  be  in  a  vertical  position  as  shown  in 
Figure  1;  i.e.,  the  deflection  of  the  cable  due  to  currents  is  small.  Any  displace¬ 
ment  of  the  cable  support  point  will  cause  the  dynamic  displacement  u(x,  t)  at  time 
t  of  an  element  of  cable  originally  located  at  x;  see  Figure  1.  Displacement 
u(x,  t)  is  of  a  form  which  satisfies  the  equation  for  the  propagation  of  waves  in  the 
cable: 


where  Pc  =  density  of  the  cable 

S  -  material  cross  section  of  the  cable 
E  =  modulus  of  elasticity  of  the  cable 

K  =  constant  of  friction  on  the  cable  due  to  the  surrounding  water 

In  an  actual  lowering  or  raising  operation,  the  length  of  the  cable,  L,  varies 
with  time.  However,  it  is  assumed  that  L  may  be  considered  constant  over  short 
periods  of  time;  i.e.,  the  net  vertical  motion  of  the  load  and  cable  does  not  influ¬ 
ence  the  dynamic  displacements  due  to  the  cable  support-point  oscillations. 

The  following  nondimensional  variables  and  parameters  may  then  be  defined: 


(4d,  e) 


The  boundary  conditions  at  the  upper  end  of  the  cable  is  the  specification  that 
u  at  x*  =  0.  The  boundary  condition  at  the  load,  when  x  =  L,  is  given  by 


M 

a 


+ 


+  2  CDpA 


du  5 u 

St"  dT 


0 


(A -2) 


where  Mq  =  dynamic  mass  of  the  array 

Cp  =  drag  coefficient  of  the  array 
A  =  horizontal  cross  section  of  the  load 
By  defining  the  parameters 


PC<L 


B  = 


CDpA 
2  M " 


Equation  A-2  may  be  reduced  to 


d^u  Bu  .  „ 

- 1  +  M  TT,  +  B 

(Bf)2  dX 


du 

dt' 


du 

St' 


=  0 


(4b,  f) 


(2) 


at  x'  =  1.0. 

The  difficulty  of  applying  this  boundary  condition.  Equation  A-2,  arises  from 
the  nonlinear  term  B  j  du/dt1  |  (du/dt'),  which  represents  the  drag  on  the  load.  To 
avoid  the  complexities  arising  from  this  nonlinear  term,  an  approximation  was  made 
in  the  reference  report  by  replacing  the  |  B u/Bt' |  term  by  (8/3tr)aJU],  where  11]  is 
the  amplitude  of  the  load  displacement,  which  is  assumed  to  be  sinusoidal.  This 
selection  results  in  the  same  energy  dissipation  when  u  is  sinusoidal  in  the  third  term 
of  Equation  2.  It  is  demonstrated  in  the  report  that  this  approximation  leads  to 
errors  on  the  order  of  20%  in  the  drag  term. 

Defining  a  normalized  displacement  amplitude  U*  equal  to  U  divided  by  |  U0[ 
and  noting  that  Uj  is  the  value  of  U1  at  the  load,  a  solution  for  U'  as  a  function  of 
x'  is  given  by 


where  y'  =  1  -  x1  and 


U1  =  U']  cosw'  y'  +  C  since' y' 


CO 


co  L 
c 


(5) 

(6b) 
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This  solution  satisfies  the  governing  equation,  Equation  1,  provided  the 
friction  of  wctcr  on  the  cable  may  be  neglected.  In  Appendix  A  of  the  reference 
report  it  is  shown  that  this  assumption  is  valid  for  the  frequency  range  of  interest. 

Substituting  Equation  5  into  Equation  2  and  incorporating  the  boundary 
conditions,  the  unknown  complex  constant  C  is  determined  as 


C  =  ±  U’jH  +  i  0  U'j ) 


and  hence,  Equation  5  reduces  to 


2 

U'  =  U\  sec<p  cos  (co1  y1  +  <p)  +  i  /3  (U\ )  tancpsinoo' 


where 


tan  <p  =  —  ,  0  *  <P  *  y 


In  requiring  that  ]  U*  at  y‘  =  1  be  equal  to  1,  U'j  is  determined  to  be 


(A-3) 


(A-4) 


(A -5) 


(y*  )2  _  cos^  (cp*  +  <p)  |  ^  +  (3*  sin^co*  s?n^2(P 

2/3^  sin^cpsin^  op'  /  cos*  (uj*  +<p) 


If  the  amplitude  of  the  dynamic  stress  is  denoted  by  I  and  a  normalized  stress 
amplitude,  I1,  is  defined  equal  to  LI/|ll0|E,  the  distribution  of  I1  is  given  by 

I‘  -  ou'1  U'i  sec<p  sin  (a,'1  y1  +  <p)  -  i  u,'/3  (Uj)^  tan<p  cos  a;’  y‘  (A-6) 

Hence  the  normalized  amplitude  of  the  maximum  dynamic  stress  j  ^ax  is  of  the  form 

(I^ox^  =  (w1)^  (Uj)^  [1  +  tan <p  (tan 'if  +  sec ^)]  \7) 

where 


2  1.222 
cos  (*'  +  <p)  |  ^  +  p  sin  oo1  sin  2o 

2 sin^<p  sin^  w*  /  cos*  (a;*  +<p) 


<P  =  arc  tan  —  ,  0  *  <p  s*  y 
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/ 


arc  tan 


tano  -  cct2<p 


Equation  7  together  with  Equations  8,  9,  and  10  was  evaluated  by  use  of  a 
computer  program  as  given  in  Appendix  C,  in  which  l^gxl  was  determined  as  a 
function  of  cc1  for  various  ranges  of  j3  and  U- 


Appendix  B 


EVALUATION  OF  THE  NORMALIZED  AMPLITUDE  OF  THE 
MAXIMUM  DYNAMIC  STRESS  AS  cu'  -  nff 


As  noted  in  the  main  text  in  the  Discussion  of  results,  the  normalized  amplitude 
of  the  maximum  dynamic  stress,  |  I^ax  |  /  is  particularly  sensitive  to  variations  of  'wc'. 
For  particular  values  of  a)1,  given  specific  values  of  j3  and/i,  it  is  desirable  to  eval¬ 
uate  the  peak  in  |^ax|  more  precisely  than  by  an  interpolation  of  the  computer 
output.  This  essentially  requires  the  derivation  of  dl^^/du;'  from  Equation  7, 
equating  this  to  zero,  and  solving  for  a;'  as  a  function  of  /9  and  fi.  Inspection  of 
Equation  7  indicates  the  complexity  of  this  derivation,  which  is  unilluminating  in 
terms  of  a  proposed  design  procedure.  As  an  alternative.  Equation  7  was  evaluated 
ir  the  limit  as  to'  approaches  ntr,  where  n  =  1,  2,  3, . . .  In  certain  cases,  this 
corresponds  to  the  peak  in  the  maximum  dynamic  stress  l^naxl*  The  analysis  was 
carried  out  in  order  to  determine  the  inaccuracies  involved  in  interpolating  the 
computer  output,  and  is  repeated  here  for  completenass. 

Given  Equations  7,  8,  9,  and  10  below,  it  is  required  to  determine  the  value 
•M^xl  as  a i*  approaches  n  ff ,  where  n  =  1,  2,  3, . . . 


(*max^  “  +  tano(ton^  +  sec'I')! 


(7) 


<U',)2  = 


coi  (u:'  +  <p) 


.2  .2  • 
zp  sin  asm  ur 


J  .  2  ,  .  20 

,  p  sin  a;  sin  lo 

T  — . . . -  ■  - 

cos^  (a;'  +<p)  J 


i  V2 


-  1  >  (8) 


a;'  IT 

Q  -  arc  tan  —  ,  0  c  *  -=r 

V  2 


(9) 


^  =  ore  tan 


y  &  (Uj)^  tone  -  cot2o 


From  Equation  9,  when  u:'  '“‘nff. 


ff  ff 

2  2 


(10) 


tone 


u;' 


nff 


(B-1) 
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cot  2 


I  -  ton 
2  tan  : 


2 


■o 


From  Equations  10,  B-l,  and  8-2, 

1 


tan  'if 


As  a.'1  -  n  n,  sin2  x'  -  0 
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1  2  .  ,..2  n  rr 

U,)  — 

2  ’  u 


;  •  •  (?) 
L’C,( 


2  2 

COS  (a.'*  +  o)  cos  o 


cos  (a:1  +  o)  "*  cos  o 


2  2 

•  *•  ^ 


By  d?fin  ng  C  (fi  sin  2c),  cos  o,  Equation  8  becomes 
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(u1,)2 


2 

cos  o 


a  ,2  .  2  .  .  2  ,  i 

24  Sin  Osm  a. 


( 1  *  C  sin2  a‘) 
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In  Equation  B-6, 
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Therefore,  from  Equation  B-6, 
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•  im 
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/  2  2  2  \ 

1  /  p  sin  2<p  cos  cp  ] 

2  4  «*2  .  2  "  J 
\  cos  cp  2p  sin  <p  / 


1  /  sin22cp  ' 


4  \  2  .2 

\cos  <psin  <p 


That  is, 


!im  (U|)^  =  1 
OJ4-  nff 


(B— 7) 


Thus,  from  Equation  B-3, 


\  i  -  (— V 

tan'!'  -  )  -5-  p2  —  -  - — 

2  M  >(*) 


But  sec 'i  -  ±y]  +  tan  ^ ,  and  since  -  ff/2  £  'I'  £  ff/2 , 


(B-8) 


sec  'i  -  +y]  +  tan  ^ 


(B-9) 


From  Equations  7,  B-l,  R-7,  and  B-9, 


2  /  _\2  I.  .  nff  .  .t.  .  .  ?.T.v^2 
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Appendix  C 


COMPUTER  PROGRAM  FOR  THE  EVALUATION  OF  THE  NORMALIZED 
AMPLITUDE  OF  THE  MAXIMUM  DYNAMIC  STRESS,  j  Z^QXj 


The  program  evaluates  the  normalized  amplitude  of  the  maximum  dynamic 
stress,  according  to  Equations  7,  8,  9,  and  10  given  in  Appendix  A  and 

is  written  for  use  on  an  IBM  1620  digital  computer.  A  flow  chart  for  the  program 
is  given  in  Figure  C-l,  followed  by  the  FORTRAN  source  program. 

Input  parameters  are  read  as  follows  in  format  (2F10.2,  2F10.8,  2F5.2,  F10.8): 


CAY 

UM 

SFR 

FMAX 

STU 

UMAX 

FZER 


4  C^A 
3»M0  k 
w  L 

=  parameter  equal  to  -rr-  =  4 

Mq 

=  increment  on  nondimensional  frequency  scale,  Awi' 

-  maximum  nondimensional  frequency,  CU^nox 

=  increment  on  input  displacement  amplitude  M*  AUo 

=  maximum  input  displacement  amplitude,  U_  I 

1  oimax 

=  initial  nondimensional  frequency,  00q 


=  parameter  equal  to 


The  computed  cjtput  is  presented  on  punched  cards  in  the  following  form.  The 
values  of  k,  p,  0,  and|  U0j  are  given  followed  by  the  input  parameters  as  defined 
above.  The  computed  values  of  J  Z^J)(  j  a-e  then  tabulated  at  each  value  of  the 
nondimensional  frequency,  according  to  format  (E15.8,  3X,  F8. 4)  fora  particular 
|  UQJ .  The  process  is  then  repeated  for  each  value  of  |  U0|  up  to  |  UQ|max,  at  which 
point  a  new  set  of  input  data  is  required. 

The  initial  frequency  is  used  as  an  input  parameter  :n  order  that  specific  ranges 
on  the  nondimensional  frequency  axis  may  be  investigated;  e.g.,  relatively  small 
increments  in  frequency  may  be  used  over  a  range  of  frequency  corresponding  to  peak 
values  in  j  There  is  a  limit  to  the  smallest  allowable  increment  in  frequency 

resulting  from  the  rounding-off  errors  inherent  to  the  program,  v*hich  results  in  | 
equaling  0  at  co'  near  to  tr.  These  errors  are  discussed  in  me  main  text  under 
Discussion  of  Results. 


Start 


Figure  C-l.  Flow  chart  for  computer  program  to  evaluate  I' 


max  l 
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n  n  n  n 


FORTRAN  SOURCE  PROGRAM 


STRESS  IN  CABLE  LOWERING  LOADS  TO  DEEP  OCEAN# 

CAY*CONSTANT  *K»  UM  =  PARAME  TER  *  SFR*STEP  IN  FREQUENCY, 

FMAX»MAX I MUM  FREQUENCY,  STU*STEP  IN  U-ZERO, 

UMAX=MAX I MUM  U-ZERO ,  FZ£R= I N I  7  I AL  FREQUENCY, 

1  READ  100, CAY, UM,SFR,FMAX,STU, UMAX, FZER 
U0*1.0 

8  BETA=CAY*UO 
PUNCH  101 

PUNCH  102, CAY, UM, BETA, UO 
FREQ*FZER 

PUNCH  111  ,CAY,UM,SFR,FMAX,STU, UMAX, FZER 
PUNCH  109 

5  PHI *A TAN ( FREQ/ UM) 

ALPH«2,*PHl 

DELT  =FREQ+PHI 

UONS* ( COS ( DEL T ) ) **2/ ( 2. *BET  A**2*S I N ( PHI )**2*SIN(FREQ)**2 ) 

TERM* ( 1 ,+ { <BETA#*2*SIN(FREQ)**2*SIN(ALPH)*»2)/COS<DtLT)**4) ) 

TERM«TERM#*0.5 

TERM*TERM-1. 

UON* UONS* TERM 
U*UON**0 , 5 

TPS  I  *  C (0,5*BETA**2,*UON)*SIN(PHI)/COS( PHI ) ) - ( COS ( AlPH ) /S IN ( ALPH ) ) 
PS  I *A  TAN ( TPS I ) 

STRS=(1,+SIN(PHI  )  /COS  (PH  I  )*<SIN(PSI  ) /COS  I  PS  I )  «-l  . /COS  <  PS  I  )  )  ) 
STRS«STRS*UON*(FREQ**2,  ) 

STR»STRS**0 • 5 
PUNCH  110, STR, FREQ 
I F ( FREQ-FMAX  )  3 ,4 , 4 

3  FREQ=FREQ+SFR 
GO  TO  5 

4  I F ( UO-UMAX ) 6 , 1  ,  1 

6  UO*UO+STU 
GO  TO  8 

100  FORMAT (2F10.2.2F1 0.8, 2F5.2.F10.8) 

101  FORMAT ( 8X ,9H  CONST ANT , 4X , 3H  UM,10X,5H  BETA,8X,7H  U-ZERO) 

102  FORMAT (3X,F10.2,3X,F10.2,3X,F10.2,3X,F10.2) 

109  FORMAT ( 7H  STRESS, 10X , 10H  FREQUENCY) 

110  FORMAT (E15.8,3X«F8.4) 

111  FORMAT (2F10.2»2F10.5»2F5.2*F10,5) 

END 
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Appendix  D 

SUMMARY  OF  DRAG  COEFFICIENTS 


INTRODUCTION 

In  an  analysis  of  the  motions  of  a  body  through  water,  whether  the  body  is 
falling  freely  or  being  lowered  by  cable,  one  of  the  most  important  effects  which 
must  be  considered  is  the  resistance,  or  drag,  experienced  by  the  body. 

The  purpose  of  this  appendix  is  to  summarize  existing  information  on  drag 
forces  and  indicate  areas  of  work  which  must  be  covered  in  order  that  such  forces 
may  be  included  in  calculating  the  motions  of  a  load  being  lowered  to  the  deep 
ocean. 


DRAG  IN  UNIFORM  FLOW 

On  the  front  of  every  solid  body  moving  through  water,  there  is  at  least  one 
point  where  there  is  no  relative  motion  between  the  water  particles  and  the  body; 
i.e. ,  there  is  a  stagnation  point.  The  pressure  at  this  point,  termed  the  dynamic 
pressure,  is  given  as 


(D-l) 


where  p  is  the  density  of  water  and  V  is  the  relative  velocity  of  the  body  to  the 
water.  It  is  convenient  to  express  the  total  drag  due  to  pressure  forces  relative  to 
this  stagnation  pressure  by  defining 


D  = 


(D-2) 


where  D  is  the  drag  force  due  to  pressure,  Cp  is  the  coefficient  of  drag,  and  S 
is  a  representative  area  of  the  body  —  either  its  frontal  or  cross-sectional  area. 
Equation  D-2  is  essentially  a  definition  of  Cp. 

The  total  drag  on  any  body  consists  of  the  "pressure  drag, "  defined  above, 
plus  drag  forces  due  to  skin  friction.  However,  for  angular  bodies  such  as  those 
envisaged  as  loads  to  be  lowered  to  the  deep  ocean  floor,  the  skin  friction  drag 
may  be  assumed  small  compared  to  the  pressure  drag. 
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Accord/  y  ;o  '•tynoids'  Similarity  Law,  the  flow  pattern  around  the  drag 
coefficients  on  two  similar  bodies  (identical  in  shape  but  dissimilar  in  size)  moving 
through  a  body  of  water  are  similar  if  their  Reynolds  numbers,  Re,  are  identical: 

R  -  ^  =  **  (D— 3) 

e  n  v 

where  V  is  the  velocity  of  the  body  relative  to  the  water,  p  is  the  density  of  the 
water,  jl  is  its  absolute  viscosity,  v  is  its  kinematic  viscosity  (ll/p),  and  d  is  a 
characterizing  dimension  of  the  body. 

Hence  it  is  possible  to  determine  the  appropriate  Cp  for  a  body  moving  through 
water  from  the  results  of  experiments  performed  on  an  identically  shaped  body  of  a 
different  scale  and  possibly  in  a  different  fluid,  provided  the  Reynolds  numbers  are 
equal. 

The  kinematic  viscosity  of  sea  water  at  normal  temperatures  and  pressures  is 
on  the  order  of  1.5  x  10“5  square  feet  per  second.  If  a  load  to  be  lowered  to  the 
deep  ocean  has  a  typical  dimension  of  15  feet,  and  moves  at  a  velocity  on  the  order 
of  1  foot  per  second,  the  Reynolds  number,  Rg,  equals  10^.  It  appears  that  relatively 
little  information  is  available  from  the  literature  on  the  variation  of  coefficients  of 
drag  at  Reynolds  numbers  greater  than  10^  to  10^.  Figure  D~1  with  inserts  show  the 
variation  of  Cp  with  Rfi  for  spheres  and  cylinders  respectively,  and  summarizes  some, 
though  by  no  means  all,  existing  data  on  the  coefficients  of  drag  applicable  to  bodies 
of  different  shapes. 

Although  objects  to  be  dropped  or  lowered  to  the  deep  ocean  floor  may  not  be 
spherical  or  cylindrical,  a  brief  investigation  of  the  dynamics  of  a  sphere  is  illu¬ 
minating.  Consider  a  body  held  stationary  in  water  and  which  is  then  allowed  to 
fall  freely.  During  the  initial  motions,  the  velocity  is  small  and  the  body  will 
accelerate  under  its  own  weight  minus  a  buoyancy  force  due  to  the  weight  of  water 
displaced,  the  drag  force  being  negligible  at  this  stage.  This  net  vertical  force  acts 
on  the  mass  of  the  body  plus  a  certain  frcction  of  its  mass  which  is  included  to 
account  for  the  water  contained  in  the  body,  if  any,  and  an  effective  mass  of  water 
to  which  accelerations  are  imparted  due  to  the  motion  of  the  body.  The  latter  terms 
are  usually  called  the  "apparent  added  mass";  the  total  mass  (body  mass  and  apparent 
added  mass)  being  termed  the  "virtual  mass."  Values  of  the  apparent  added  mass  vary 
from  40%  to  150%  of  the  mass  of  the  body. 

As  the  velocity  increases  from  zero,  the  drag  force  opposing  the  motion  becomes 
significant,  and  at  a  particular  time  t  =  tj  this  force  is  given  by 

fd  *  co(ipVi2)s  <Dwt) 


70 


where  Vj  is  the  velocity  at  t  -  tj ,  and  Cp  is  the  appropriate  coefficient  of  drag. 
After  a  given  time  the  velocity  of  the  falling  body  attains  a  terminal  velocity,  Vy , 
in  which  condition  the  drag  force  balances  the  body  weight  minus  the  buoyancy  force; 
i.  e. , 

<wb-fb>  -  fd  =  <cd)t(tpVt2)s  (d-5) 


where  (C^)^  is  the  coefficient  of  drag  at  the  Reynolds  number  conesponding  to  a 
velocity  of  Vy.  Equation  D-5  may  be  rewritten  as 


(WB  -  FB>  ’  FD 


(D-6) 


The  function  f(VjL/ju)  is  not  known  and  cannot  be  defined  analytically,  and 
Equation  D-6  cannot  be  solved  explicitly  for  the  terminal  velocity,  Vy,  without 
the  prior  assumption  of  a  particular  Cq. 

However,  starting  from  zero  initial  velocity,  it  is  possible  to  determine  the 
motion  of  a  particular  body  by  considering  the  acceleration  and  velocities  attained 
over  small  increments  of  time.  A  simple  computer  program  was  written  to  accomplish 
this.  At  t  =  0  the  velocity  is  zero,  there  is  no  drag  force,  and  the  body  will  accel¬ 
erate  under  the  force  (Wg-Fg).  At  t  =  t]  the  velocity  is  finite  and  the  appropriate 
Reynolds  number  may  be  calculated  together  with  the  corresponding  Cq  .  For  the 
purpose  of  these  calculations,  Cq  was  specified  at  increments  of  Rfi,  and  a  simple 
interpolation  was  made  to  determine  the  specific  Cp  corresponding  to  Re  at  t  =  tj  . 
Hence  the  drag  and  out-of-balance  force  may  be  calculated  at  t  =  tj  together  with 
the  instantaneous  acceleration  at  this  point.  The  process  can  be  repeated  to  deter¬ 
mine  the  velocity  of  the  body  at  time  increments  from  t  =  0  to  t  =  T,  where  T  is  the 
time  taken  to  attain  a  terminal  velocity. 

For  an  up-to-date  complete  treatment  of  hydrodynamic  drag,  the  excellent 
treatise  prepared  and  published  by  Dr.  Sighard  F.  Hoerner^  should  be  consulted. 
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Reynold*  Numbtr,  R#  -p— 

Voriation  of  Cq  with  Reynolds  number 


Appendix  E 

SUMMARY  OF  ADDED  MASS  COEFFICIENTS 


The  concept  of  added  mass  is  well  known  in  fluid  mechanics.  The  physical 
explanation  of  this  phenomenon  is  that  when  a  body  is  subjected  to  an  unbalanced 
force,  not  only  must  the  mass  of  the  body  be  accelerated,  but  also  that  of  the 
added  fluid  mass  surrounding  the  body.  The  ratio  of  this  added  fluid  mass  to  the 
body  mass  is  i! _ ;I  J:^  mas;  coefficient,  Cm. 

The  added  mass  depends  on  the  dimensions  and  shape  of  the  body  and  the 
density  and  viscosity  of  the  fluid.  In  general,  measurements  of  the  apparent  added 
mass  have  been  obtained  under  two  fundamentally  different  flow  situations.  In  one 
the  motion  is  oscillatory  in  that  an  immersed  body  is  vibrated.  In  the  other  the 
motion  is  unidirectional  in  that  an  immersed  body  is  accelerated  recti  linearly.  The 
exact  analytical  description  of  fluid  resistance  to  the  acceleration  of  an  arbitrarily 
shaped  submerged  solid  is  not  known,  hence  the  exact  added  mass  coefficient  for 
various  shapes  of  objects  is  not  known.  The  following  reports  are  the  results  of 
different  experiments  under  different  conditions,  but  they  are  quite  consistent: 

T.  E.  Ste Ison  and  F.  T.  Mavis,  ^  E.  Silberman,  ?  T.  Sarpkaya,  ®  N.  L.  Ackermann 
and  A.  Arbhabhirama,  ^  O.  C.  Zienkiowicz  and  B.  Nath.  *0  A  summary  of  most  of 
the  important  results  obtained  from  these  references  is  presented  in  Figure  E-l, 

The  coefficients,  have  been  arranged  in  terms  of  a  common  dimension, 
namely  the  ratio  of  the  added  mass  to  the  mass  of  fluid  displaced.  The  results  obtained 
from  oscillatory  motion  are  as  follows: 

1.  Spheres:  Cm  -  0.51.  This  compares  with  a  value  of  0. 50  obtained  from 
ideal  fluid  theory  for  rectilinear  motion. 

2.  Cubes:  <!„,  =  0.67  ("broodside-on"  or  "edge-on"). 

3.  Circular  Cy linden:  See  Figure  E-l.  The  abscissa  is  the  ratio  of  length  to 
diameter.  The  motion  is  in  the  direction  perpendicular  to  the  circular 
cross  section. 

4.  Rectangular  Plates:  See  Figure  E-l.  The  abscissa  is  the  ratio  of  length  to 
width.  The  motion  is  in  the  broadside-on  direction.  The  ratio  of  thickness 
to  width  is  limited  to  values  less  than  0.04. 

5.  Square  Prisms:  See  Figure  E-l.  The  abscissa  is  the  ratio  of  length  to  width 
of  the  square  sides.  The  motion  is  in  the  direction  perpendicular  to  the 
square  cross  section. 
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6.  Symmetrical  Lenses:  See  Figure  E-l.  These  lenses  are  two  intersecting  or 
separated  spheres: 


— « - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 
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The  abscissa  is  the  ratio  of  B/R.  The  motion  is  in  the  direction  shown. 

7.  Two  Parallel  Rectangular  and  Square  Plates:  See  Figure  E-l.  The  abscissa 
is  the  ratio  of  spacing  between  two  plates  to  width,  where  the  spacing  is 
measured  from  center  to  center  of  the  plates.  The  ratio  of  length  to  width 
of  the  plates  ara  over  17  to  1.  The  motion  is  in  a  direction  parallel  to  the 
thickness  of  the  plates. 

For  all  of  the  oscillatory  motion  cases  the  experiments  were  conducted  at  low 
velocities  and  high  accelerations.  Thus  the  total  resistance  force  to  the  moving 
object  is  largely  due  to  the  added  mass  which  is  dependent  on  the  acceleration.  At 
higher  velocities,  the  total  resistance  force  is  due  to  a  velocity -dependent  drag  term 
as  well  as  to  the  added  mass  term.  That  port  of  the  resistance  to  motion  due  to  viscous 
and  form  drag  and  that  part  due  to  added  mass  are  difficult  to  separate.  Sielson  and 
Movis^  and  Silberman?  realized  the  difficulties.  From  experiments  on  a  sphere  the 
measured  added  mass  increased  by  approximately  1%  above  the  values  obtained  from 
ideal  fluid  theory.  Thus  it  was  concluded  that  viscosity  did  not  seriously  affect  the 
experimental  values  for  the  added  mass. 

A  recent  method  (Zienkiewicz  and  Noth  10)  of  measuring  the  added  moss  is 
worth  mentioning  here,  thing  an  electric  analogy  method,  the  virtual  mass  as  well 
as  the  pressure  distribution  around  a  rigid  body  accelerating  in  an  incompressible 
fluid  can  be  determined.  In  the  following  table,  the  results  are  compared  with  the 
known  odded  moss  coefficients  obtained  from  other  sources.  Agreement  is  excellent. 


Added  Moss  G>efficient 

Object 

Obtained  by 

Zienkiewicz  10 

From  Indicated  Source 

Infinitely  long  vertical  plate 

1.03 

1.04  (Riabouchintky  1®) 

Infinitely  long  cylinder 

0.98 

1.00  (H.  Iambi  1) 

Thin  circular  disc 

0.61 

0.636  (H.  Lombn) 

Gibe 

0.62 

0.67  ($telson6) 

.  . .  . * 
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These  values  are  measured  for  an  infinitely  large  submergence  depth.  At  small 
depths  the  measured  added  mass  decreased.  This  agrees  with  the  physical  explanation 
of  the  added  mass  phenomenon.  The  experiments  were  conducted  for  small -amplitude 
motions  and  no  separation  occurred,  hence  the  boundary  effect  is  not  considered. 

This  method  can  be  used  for  rotary  acceleration  of  o  body.  This  experimental  method 
has  important  implications  since  it  can  be  set  up  in  a  deep  ocean  simulating  tank  to 
measure  the  added  mass  as  well  as  the  pressure  distribution  of  any  arbitrarily  shaped 
object  under  translational  or  rotary  motion. 

For  the  unidirectional  motion,  the  viscous  and  boundary  effects  must  be 
considered.  Experiments  for  this  type  of  motion  have  been  conducted  for  a  number 
of  bodies,  but  only  that  for  spheres  will  be  cited.  Arbhabhirama^  found  that  when 
the  ratio  of  the  diameter  of  a  sphere  to  the  diameter  of  a  fluid  filling  a  concentric 
spherical  shell  is  0.259,  which  is  similar  to  a  sphere  oscillating  in  an  infinite  fluid, 
the  added  mass  is  found  to  be  1.03  times  the  added  mass  obtained  from  potential  flow. 

In  summary,  although  some  data  is  available  on  added  mass  coefficients  in 
oscillatory  flow,  most  of  the  experiments  hove  been  conducted  at  small  scale  and 
within  the  low  Reynolds  number  regime.  As  an  example,  the  following  estimate  of 
the  added  mass  of  a  complicated  frame  structure  such  as  the  STU  described  below,  is 
cited.  Theoretically  the  oscillatory  motion  of  a  load  being  lowered  to  the  ocean 
floor  and  suspended  by  a  cable  is  a  damped  simple  harmonic  motion.  If  the  cable  is 
considered  to  be  elastic,  the  equation  of  motion  is 

Mx  +  Cx  +  kx  =  0 
a 

where  M  is  the  virtual  mass  of  the  load,  x  is  the  elongation  of  the  cable,  k  is  the 
ratio  of  the  restraining  force  to  the  elongation  of  the  coble,  and  C  is  the  coefficient 
of  damping.  The  solution  of  this  equation  is  quite  complicated,  but  the  period  of 
oscillation  is  the  same  as  for  simple  undamped  harmonic  motion: 

Mq  x  +  k  x  r  0 
Hence,  the  period  T  =  2tt(MQA). 

On  13  April  1965,  o  Submersible  Test  Unit  (STU)  was  lowered  by  this  Laboratory 
to  the  ocean  floor  to  o  depth  of  2,500  feet  using  a  1.3-inch-diameter  polypropylene 
cable.  The  cable  tensions  were  recorded  os  a  function  of  time  from  the  start  o? 
lowering  operation.  The  curve  of  the  graph  (Figure  35)  shows  the  oscillatory  motion 
of  the  STU,  consequently  the  overage  period  of  9.8  seconds  wos  obtained  while  the 
average  tension  of  the  coble  in  this  interval  (between  10  and  12  minutes  as  marked  on 
the  figure)  was  3,400  pounds.  The  breoking  strength  of  the  cable  is  45,000  pounds 
C Braided  Rope  and  Cordage  Catalog,  ”  Somson  Cordoge  Works,  Boston,  Mass. ).  From 
the  percent  load  of  breaking  strength  versus  the  percent  elongation  curve  of  polypro¬ 
pylene  cobles,  the  corresponding  percent  elongation  of  4.0  is  obtained.  Since  the 
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average  cable  length  in  this  interval  is  586  feet,  the  elongation  of  the  cable  is 
0  04  x  586  =  23.4  feet.  The  restoring -force  constant  was  determined  to  be  145  pounds 
per  foot  3,400/23.4).  The  virtual  mass  of  the  coble  assembly  was  found  to  be 

1 1,350  pounds;  i.  e. ,  Mq  -  (T^/4^)g  =  [  (9.8)2  (145)/42  j  32.2.  Of  this  1 1,350  pounds, 
950  pounds  is  due  to  the  static  weight  of  the  suspended  cable.  The  net  weight  of  the 
STU  in  wcAer  is  5,500  pounds.  Thus  the  added  mass  coefficient,  C^,  is  determined 
to  be  slightly  less  than  2.0  (i.e.,  10,700/5,500). 
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(  See  text  listing  in  Appendix  E  for  abscissa  units) 

Figure  E-l.  Added  mass  coefficient  as  a  function  of  circular  frequency 
for  constant  geometric  ratios. 


REFERENCES 


1.  Arthur  D.  Little,  Inc.  Technical  Report  no.  1370863  on  Project  Trident,  Bureau 
of  Ships  Contract  NObsr-81564:  Stress  analysis  of  ship-suspended  heavily  loaded 
cables  for  deep  underwater  emplacements.  Cambridge,  Mass.,  Aug.  1963. 
AD-418028. 

2.  Marine  Advisers,  Inc,  Report  on  Contract  NBy-32206:  The  motions  of  a  moored 
construction-type  barge  in  irregular  waves  and  their  influence  on  construction 
operation,  by  P,  Kaplan  and  R.  R.  Putz.  La  Jolla,  Calif.,  Aug.  1962. 

3.  U.  S.  Naval  Civil  Engineering  Laboratory.  Technical  Note  N-604:  The 
engineering  applications  of  a  report  entitled  "The  motions  of  a  moored  construction- 
type  barge  in  irreguiar  waves  and  their  influence  on  construction  operation,"  by 

W.  J.  Pierson  and  P.  Holmes.  Port  Hueneme,  Calif.,  Aug.  1964. 

4.  U.  S.  Navy,  David  Taylor  Modei  Basin.  Report  1221:  Theoretical  analysis  of 
the  effect  of  ship  motion  on  mooring  cables  in  deep  water,  by  L.  F.  Whicker. 
Washington,  D.  C. ,  Mar.  1958. 

5.  S.  F.  Hoerner.  Fluid-dynamic  drag,  2nd  ed.  Author,  148  Bustead  Drive, 

Midland  Park,  N.  J.,  1965. 

6.  T.  E.  Stetson  and  F.  T.  Mavis.  "Virtual  mass  and  acceleration  in  fluids," 
American  Society  of  Civil  Engineers,  Transactions,  vol.  122,  1957,  paper  no.  2870, 
pp.  518-525. 

7.  E.  Silberman.  Discussion  of  the  paper,  "Virtual  mass  and  acceleration  in  fluids," 
by  T.  E.  Stetson  and  F,  T.  Mavis.  American  Society  of  Civil  Engineers,  Transactions, 
vol.  122,  1957,  pp.  527-529. 

8.  T.  Sarpkaya.  "Added  mass  of  lenses  and  parallel  plates,"  American  Society  of 
Civil  Engineers,  Proceedings,  Journal  of  the  Engineering  Mechanics  Division, 
vol.  86,  EM3,  June  1960,  pp.  141-152, 

9.  N.  L.  Ackermann  and  A.  Arbhabhirama.  "Viscous  and  boundary  effects  on 
virtual  maw,"  American  Society  of  Civil  Engineers,  Proceedings,  Journal  of  the 
Engineering  Mechanics  Division,  vol.  90,  EM4,  Aug.  1964,  pp.  123-130. 

10.  O.  C.  Zienkiewicz  and  B.  Nath.  "Analogue  procedure  for  determination  of 
virtual  mass,"  American  Society  of  Civil  Engineers,  Proceedings,  Journal  of  the 
Hydraulics  Division,  vol.  90,  HY5,  Sept.  1964,  pp.  69-81, 

11.  H.  Lamb.  Hydrodynamics,  6th  ed.,  1932.  London,  Cambridge  University 
Press;  New  York,  Dover  Publications. 


MATHEMATICAL  NOTATIONS 


A  Cross-sectional  area  of  the  load  in  the  direction  of  motion 
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p  Density  of  seawater 
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MATHEMATICAL  NOTATIONS 
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